On the nature of radio pulsars with long periods 
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ABSTRACT 

It is shown that the drift waves near the light cylinder can cause the modu- 
lation of emission with periods of order several seconds. These periods explain 
the intervals between successive pulses observed in "magnetars" and radio pul- 
sars with long periods. The model under consideration gives the possibility to 
calculate real rotation periods of host neutron stars. They are less than 1 sec for 
the investigated objects. The magnetic fields at the surface of the neutron star 
are of order 10 11 - 10 13 G and equal to the fields usual for known radio pulsars. 

Subject headings: pulsars: individual (PSR J2144-3933, PSR J1847-0130, PSR 
J1814-1744)-stars: magnetic fields 



1. Introduction 

According to the plasma model pulsar radio emission is generated in the electron- 
positron plasma, which appears by the avalanche process of (e + e~) pair production. The 
accomplishment of this process requires fulfilment of the following condition: 



Bschw > B s > B d i 



(1) 
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Here B s = 3.2 x 10 19 \fp~P is pulsar surface magnetic field inferred from observational 
values of spin period and period derivative, B Sc hw = m 2 e c' feh 4.4 x 10 13 G is Schwinger 
limit i.e. the magnetic field at which the electron cyclotron energy equals the electron 
rest mass energy and B^i is the magnetic field corresponding to the so-called "death line" 
(see Fig. 1) (Young et al. 1999). The various configurations of the surface magnetic field 
correspond to different death lines. They depend not only on magnetic field configuration, 
but significantly on whether the origin of gamma quanta, which are responsible for pair 
production, is curvature radiation or inverse Compton scattering (Zhang & Harding 2000). 
As a mechanism of producing of gamma quanta we took curvature radiation with sunspot 
configuration of magnetic field. In this case "death line" is defined by the condition that the 
potential drop across the gap required to produce enough pairs per primary to screen out the 
parallel electric field is larger than the maximum potential drop available from the pulsar. 
The quantity B d i can be solved from the following equation (Chen & Ruderman 1993): 

71og5 ( «-131ogP = 78 (2) 

Recently there were discovered three long period radio pulsars PSR J2144-3933 (Young 
et al. 1999), PSR J1847-0130 (McLaughlin et al. 2003) and PSR J1814-1744 (Camilo et 
al. 2000) (see Table 1), that break the condition above (the first one breaks the right-hand 
side of the inequality and the other two break the left). PSR J2144-3933 is distinguished 
by some other characteristics. It has the lowest spindown luminosity (E = 4ir 2 IP/P 3 m 
3.2 x 10 28 erg • s _1 ) of any known pulsar. The beaming fraction (that is, the fraction of 
the celestial sphere swept across by the beam) is also the smallest, W « 1/300. On the 
other hand, PSR J1847-0130 (McLaughlin et al. 2003) and PSR J1814-1744 (Camilo et al. 
2000) are isolated radio pulsars having the largest, " magnet ar" -like, inferred surface dipole 
magnetic fields yet seen in the population: 9.4 x 10 13 G and 5.5 x 10 13 G, respectively. These 
pulsars show apparently normal radio emission in the regime of magnetic field strength 
(B s > Bschw), where plasma models predict no emission should occur. However, the nature 
of the Schwinger limit is not clear and is the subject of long term discussions. Baring & 
Harding (2001) proposed that for extremely strong fields, photon splitting dominates over 
pair production and the corresponding death line is a function of both field strength and 
period. This is definitely neither a certain nor a hard limit. Furthermore, Usov (2002) has 
argued that since only one photon polarization mode splits, the onset of photon splitting 
does not suppress the production of pairs. It is very important to investigate this interesting 
problem, but it is beyond the framework of our model, especially since we argue that magnetic 
fields of neutron stars in our investigation are less than the Schwinger limit. 

A model that explains the phenomenon of radio emission from these pulsars and all 
their anomalous properties does not exist. 
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An important feature of our model, which provides a natural explanation of most of 
the properties of these pulsars, is the presence of very low frequency, nearly transverse drift 
waves propagating across the magnetic field and encircling the open field line region of the 
pulsar magnetosphere (Kazbegi et al. 1991b, 1996). These waves only periodically change 
the direction of the radio emission and are not directly observable. We give a description of 
this model in Sections 2, 3 and 4. Some estimates of the rotational and angular parameters 
of these pulsars are given in Section 5. We discuss the obtained results in Sections 6 and 7. 

2. Emission model 

As mentioned the pulsar magnetosphere is filled by a dense relativistic electron-positron 
plasma. The (e + e~) pairs are generated as a consequence of the avalanche process [first 
described by Sturrock (1971)], and they flow along the open magnetic field lines. The 
plasma is multicomponent, with a one-dimensional distribution function (see Arons 1981a, 
Fig.l), containing the following: 

(i) electrons and positrons of the bulk of plasma with a mean Lorentz factor of 7 P and 
density n p ; 

(ii) particles of the high-energy 'tail' of the distribution function with 7 t and n t , stretched 
in the direction of positive momentum; 

(iii) the ultrarelativistic (7;, ~ 10 6 ) beam of primary particles with so-called 'Goldreich- 
Julian' density (Goldreich & Julian 1969) n b « 7 x 10- 2 B S P- 1 (R / 'R) 3 [cm- 3 ] (where P is 
the pulsar period, Rq ~ 10 6 cm is the neutron star radius, B s is the magnetic field value at 
the stellar surface and R is the distance from the center of the neutron star). This density 
is much less than the density of secondary pairs n p . 

Let us note that such parameters are only possible under the assumption that strong 
and curved multipolar fields persist in the pair-production region, near the stellar surface 
(Arons 1981b; Machabeli & Usov 1989). 

Such a distribution function should generate various wave-modes under certain condi- 
tions. The waves then propagate in the pair plasma of the pulsar magnetosphere, trans- 
form into vacuum electromagnetic waves as the plasma density drops, enter the interstellar 
medium and reach an observer as pulsar radio emission. This waves leave the magnetosphere 
propagating at relatively small angles to the pulsar magnetic field (Kazbegi et al. 1991a). 

Particles moving along the curved magnetic field undergo drift motion transversely to 
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the plane where the field line lies. The drift velocity can be written as 

u = c^v^/ujbRc (3) 

where ujb = eB/mc, B = B s (Rq/R) 3 . R c is the curvature radius of the dipolar magnetic 
field line, 7 is the relativistic Lorentz factor of a particle and is the particle velocity along 
the magnetic field. Here and below the cylindrical coordinate system (x, r, (p) is chosen, with 
the x-axis directed transversely to the plane of the field line, while r and (p are the radial 
and azimuthal coordinates, respectively. 

Generation of radio emission is possible if at least one of the following resonance condi- 
tions is fulfilled: 

LOB 

cyclotron instability u — k v v v — k x u = (4) 

Cherenkov instability u — k^v^ — k x u = (5) 

These conditions are very sensitive to the parameters of the magnetospheric plasma, particu- 
larly to the value of the drift velocity (see eq. [3]), and hence to the curvature of the magnetic 
field lines (Kazbegi et al. 1996). 

It should be noted that, in the absence of drift motion, the ordinary Cherenkov interac- 
tion implies that the phase velocity of a wave equals the velocity of the particles in both value 
and direction. In other words, an observer moving with the same velocity as the particles, 
should detect the same phase of the wave for a sufficiently long time (r ~ 1/T ). This is, 
however, impossible for a wave propagating transversely to the particles' velocity. On the 
other hand, the drift velocity (equation (3)) is directed along the phase velocity of such a 
wave. This allows wave-particle resonant interaction. 



3. Generation of drift waves 

It has been shown (Kazbegi et al. 1991b, 1996) that, in addition to the waves mentioned 
above (the characteristic frequencies of which fall into the radio band) propagating with small 
angles to the magnetic field lines, very low frequency, nearly transverse drift waves can be 
excited. They propagate across the magnetic field, so that the angle between k and B is 
close to 7r/2. In other words, k±/k^ ^> 1, where k± = (k 2 + k^) 1 / 2 . Assuming ^{uj/ujb) <C 1, 
(u a /c) 2 <^ 1, k^/kx <^ 1 and k r — > 0, we can write the general dispersion equation of the 
drift waves in the following form (Kazbegi et al. 1991a,b, 1996): 



^ to J v^t 



ecu — k^Vp — k x u a dp uj 
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^ uj J uj - k^Vp - k x u a dp to 2 



CO 



2 



+ J2^f "° /C t B -^iv\ =0 (6) 
^ w J w - fc^ - k x u a dp I 



where a denotes the sort of particles (electrons or positrons) and io 2 a = 47rn a e 2 /m, f a is the 
distribution function and p is the momentum of the plasma particles. 

Let us assume that 

oj = k^Vp + k x u h + iT (7) 

where is the drift velocity of the beam particles (see eq. [3]). In the approximation k^v^ <C 
k x u b and k 2 x <C 0^/7^ c 2 , the imaginary part can be written as 

r = Imuj » (^j V2 (Jlj ' k x u b (8) 

According to equation (7), the frequency of a drift wave can be written as 

co>o = Rgoj = k v v v + k x Ub ~ k x Ub (9) 

Drift waves propagate across the magnetic field and encircle the region of the open field 
lines of the pulsar magnetosphere. They draw energy from the longitudinal motion of the 
beam particles, as in the case of the ordinary Cherenkov wave-particle interaction. However, 
they are excited only if k x Ub 7^ 0, i.e., in the presence of drift motion of the beam particles. 
Note that these low-frequency waves are nearly transverse, with the electric vector being 
directed almost along the local magnetic field. Let us note that although k^v^ <^ k x Ub for 
the drift waves, there still exists a nonzero k v . It appears that growth rate (equation (8)) is 
rather small. However, the drift waves propagate nearly transversely to the magnetic field, 
encircling the magnetosphere, and stay in the resonance region for a substantial period of 
time. Although the particles give a small fraction of their energy to the waves and then 
leave the interaction region, they are continuously replaced by the new particles entering 
this region. The waves leave the resonance region considerably slower than the particles. 
Hence, there is no sufficient time for the inverse action of the waves on the particles. The 
accumulation of energy in the waves occurs without quasi-linear saturation. The amplitude 
of the waves grows until the nonlinear processes redistribute the energy over the spectrum. 
As was demonstrated by Kazbegi et al. (1991b), the strongest nonlinear process in this case 
is the induced scattering of waves on plasma particles. Therefore, the growth of the drift- 
wave amplitude continues until the decrement of the nonlinear waves T^l becomes equal to 
the linear decrement T. As a result, one obtains quasi-regular configurations of drift waves. 
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Generally, the nonlinear scattering pumps the wave energy into the long-wavelength domain 
of the spectrum. 

cP 

Amax ~ T LC = — (10) 

Here r LC is the radius of the light cylinder. 

According to equations (9), (10) and (3), the period of the drift waves can be written 

as: 

2vr 2vr A e BP 2 . , 

to dr k x u h u b 4n 2 mc j b 
It appears that the period of the drift wave is of the order of several seconds. It is possible 
to determine the relationship between P^ r , the derivative and the rate of slowing down of 
the neutron star from equation (11) 

eB 

P dr = l D P p (12) 
27r 2 mc7f, 

For the considered values of the parameters we obtain Pd r ~ 10P This relation is kept during 
the entire life of the pulsar, until it stops emitting. 



4. Mechanism of field line curvature change 

Let us assume that a drift wave with the dispersion defined by equation (7) is excited 
at some place in the pulsar magnetosphere. It follows from the Maxwell equations that 
B r = E v (k x c/u), hence B r ^> for such a wave. Therefore, excitation of a drift wave 
causes particular growth of the r-component of the local magnetic field. 

The field line curvature p c = 1/R C is defined in a Cartesian frame of coordinates (x, y, z) 
(where z axis is directed perpendicular to the plane of the field line) as 

-3/2 

il-n 

(13) 



Pc = 



ax 



dx 2 



where dy/dx = B y /B x . Using (V • B) = and rewriting equation (13) in cylindrical coordi- 
nates, we obtain 

1P,_11_P|9P, 
Pc rB rBB 2 dp 1 J 



Here B = (B 
equation (14) 



1/2 _ "l 

+ B 2 ) ps B v [l + B 2 /2B 2 ]. Assuming that k^r > 1, we obtain from 



B r 

Dy 



(15) 
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From equation (15) it is clear that even a small change of B r causes a significant change of 
p c . The variation of the field line curvature can be estimated as 



Kr-FT ( 16 ) 



It follows that even a drift wave with a modest amplitude B r ~ AB r ~ O.OIB^ alters the 
field line curvature substantially, Ap c / p c ~ 0.1 

Since radio waves propagates along the local magnetic field lines, curvature variation 
causes a change of the emission direction. 



5. The model 

There is unequivocal correspondence between the observable intensity and a (the angle 
between the observer's line of sight and the emission direction (see Fig. 2)). The maximum 
of the intensity corresponds to the minimum of a. The period of a pulsar is the time interval 
between the neighboring maximums of the observable intensity (minimums of a). According 
to this fact, we can say that the observable period is representative of the value of a and as 
it appears below, it might differ from the spin period of the pulsar. 



cosa = A-K (17) 

where A and K are unit guide vectors of observer's and emission axes, respectively. In the 
spherical coordinate system (r,tp,0), combined with the plane of the pulsar rotation, these 
vectors can be expressed as: 

A = (1,0, 5) (18) 

K = (i,nt,p) (19) 

where Q = 2-7T jP is the angular velocity of the pulsar. 5 is the angle between the rotation 
and observer's axes, and (3 is the angle between the rotation and emission axes (see Fig. 2). 

From equations (17), (18) and (19), it follows that: 

a = arccos (sin 5 sin (3 cos fit + cos S cos (3) (20) 

In the absence of the drift wave (3 = fi = const, and consequently the period of a is equal 
to 2ir/tt. 

According to equation (16), in the presence of the drift wave, the fractional variation 
Ap c / p c is proportional to the magnetic field of the wave B r , which is periodically changing. 
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So j3 = (3 (t) is harmonically oscillating about (3q with an amplitude A/3 = Ap c /p c and rate 
Udr = ■ Thus, we can write that 

P = Po + APsm(u dr t + <p) (21) 

According to equations (20) and (21,) we obtain 

a = arccos (sin S sin (/3 + A/3 sin (c^ r t + </?)) cos fit + cos S cos (/3 + A/3 sin {u) dr t + </?))) 

(22) 

« min = | (A - 6) + A/? sin (2ttA;^ + V ) | (23) 

a min is the minimum of a after k revolutions of the pulsar 1 . The parameters of the pulse 
profile (e.g. width, height etc.) significantly depend on what the minimal angle would be 
between the emission axis and the observer's axis while first one passes the other (during 
one revolution). If the emission cone does not cross the observer's line of sight entirely (i.e., 
the minimal angle between them is more than cone angle see inequality (24a)), then we 
cannot observe the pulsar emission. On the other hand, inequality (24b) defines condition 
that is necessary for emission detection. 

«min > $ (24a) 

oU < * (24b) 

Hence, for some values of parameters Q, u^ r , f3, A/3, 5, (p and d (Set A) it is possible to 
accomplish the following regime: after every k = m turn, the minimal value of a (a^ in ) 
satisfies condition (24b) while for intervening values oi k (1 < k < (m — 1), where k and m 
are positive integers), a™ in satisfies condition (24a). In that case the observable period P b s 
does not represent the real pulsar spin period, but is divisible by it. 

P obs = mP (25) 
Pobs = mP, (26) 



Hence, 



where P is the pulsar spin period. 

The dipolar magnetic field strength on the neutron star surface can be written as: 

B = 3.2 x 10 19 \/AP (27) 



The detection moment of any pulse is taken as the zero point of the time reckoning 
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From equations (25), (26) and (27), it follows that 

B = — (28) 
m 

After inserting equations (28) and (25) in equation (2), we obtain: 

7 log B - 13 log P = (7 log B obs - 13 log P obs ) + 6 log m > 78 (29) 

Then 

6 log m > 78 - 7 log B obs + 13 log P obs (30) 

It can be verified that there exists value for m that satisfies equation (30) and the following 
condition simultaneously 

B b 

B = < Bschw (31) 

m 

Thus, it is possible to fulfil the conditions necessary for (e + e~) pair production for some 
values of m. 

There is a common problem for pulsars II and III presented in Table 1 (their surface 
magnetic field strength exceeds the quantum critical value B > Bs c hw), whereas for PSR 
J2144-3933 the problem appears in a different way. Equation (29) is not fulfilled. After 
substituting B obs and P obs from Table 1 in equation (30), we obtain: 

61ogm>4 (32) 

Thus, if equation (32) is accomplished, it is possible for PSR J2144-3933 to generate radio 
emission. 

For better estimation of m we can use observational data for beaming fractions. From 
Figure 2 appears that the pulse width can be expressed as: 

W = P x 2sin#/27rsin(5, (33) 

after inserting equation (25) in equation (33), we obtain: 

mW/P obs = sin #/tt sin 5 (34) 

As was mentioned above, to accomplish the described regime (eq.[25]), (the angle between 
the observer's line of sight and the emission direction after one revolution from that moment 
when they were coincident (a = 0, Fig. 3)) a™ in must exceed 9. Since 

a l mm = A(3(2nP/P dr ), (35) 
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if we assume that (3q = 5, then we get Pdr = 2P &s = 2mP and 



0, 



max 



A/3(sin7r/m); 



(36) 



if we substitute this equation in equation (34), we obtain: 



W/P bs = A / 3(sin7r/m)/m7rsin(5 



(37) 



Here the left-hand side is known from the observations. Equation (37) gives us the ability 
to estimate the angular parameters of pulsars for the given values of m. 

If we consider these pulsars in the framework of our model, their parameters (e.g. spin 
and angular parameters) will get new 'real' values, shown in Table 2. 

According to the obtained results, the considered pulsars are placed on (P-B) diagram 
as shown in Fig. 4 

Thus, we developed the theoretical model of pulsar emission, in the framework of which 
we explained all specific features of pulsars presented in Table 1. 



It should be noted that this model is applicable to the entire population of pulsars, but 
the effects caused by drift waves are different depending on the values of the parameters (Set 
A). In the case of large A/3 (A[3 > ■$), the most interesting effect is the lengthening of the 
observable period (see eq.[25]) which is accomplished only when Pd r is divisible by P to high 
accuracy. It explains the lack of such kinds of pulsars. 

In the case of small A/5 (A/3 < the observable period does not increase (except 
for \j3o — 5\ ~ A/3 ), but some other interesting effects appear, such as drifting subpulses 
(Kazbegi et al. 1991c) and period and period derivative oscillation phenomenon, which is 
observed in PSR B1828-11 (Stairs et al. 2000) and PSR B1642-03 (Shabanova et al. 2001). 
Some authors (Jones & Andersson 2001; Link & Epstein 2001; Rezania 2003) have proposed 
different models to explain this phenomenon within the framework of free precession of the 
neutron star. As shown by Shaham (1977) and Sedrakian et al. (1999) the existence of 
precession in the neutron star is in strong conflict with the superfluid models for the neutron 
star interior structure. Therefore, we can declare that there does not yet exist a self-consistent 
explanation of this fact. We plan to study this problem in detail in a forthcoming paper. 

If Pdr is not divisible by P, then the observed intensity must be modulated with the 
period of the drift wave. It is impossible to get such variations with integrated pulse profiles. 



6. Discussion 



- 11 - 



Deviations of integrable pulse intensities damp each other. The only possible way to prove 
this scenario is by single pulse observations. Such observations really show intensity varia- 
tions (Karastergiou et al. 2001). Although they do not have have a harmonic nature (this 
is due to various noises and insufficient resolution), it benefits our model. So if experiments 
are modified to evolve the oscillating component, it would help validate our theory. 

Let us consider pulsars with very short periods. As mentioned above, drift waves arise 
in the vicinity of the light cylinder. The shorter the pulsar spin period, the smaller is the 
radius of the light cylinder, and consequently, the larger is the magnetic field value in the 
wave generation region [B^ m B — B s (Rq/ R) 3 ). Thus, if we take this consideration into 
account, from equation (16) it follows that for pulsars whose period is much less than 0.1s, 
the amplitude of the oscillation in the emission direction would be so small (A/3 < 1°), that 
the presence of drift waves would not cause any significant effect. 

We believe that if the scenario defined above (enlarging of the observable period due to 
emission direction changing by drift waves) is accomplished then there must be simultaneous 
modulation of the radio and X-ray emission. Some evidence of this fact is detected radio 
emission from SGR 1900+14 (Shitov 1999) and AXP 1E2259+586 (Malofeev & Malov 
2001). Higher radio frequency emission from these magnetars has not been observed yet, 
but discussion about this issue differs from objectives of our paper. However, detection 
of both types of emission is rare event, because they are generated on different altitudes, 
radio and X-ray emission propagate in different directions. This implies that the pulsars 
considered in this paper do not show X-ray emission. 



7. Conclusions 

After these considerations we can divide radio pulsars into the following groups, which 
are listed with their main requirements: 

1. Rapidly rotating pulsars, for which A/3 is too small. None of the mentioned effects 
there should exist for these. 

2. Pulsars with A/3 < $ and (P dr — P)/P dr « 1. In this case, period, period derivative 
and pulse shape oscillation should appear. In the case of low accuracy of equality between 
P dr and P, subpulse drift can be observed (Kazbegi et al. 1991a; Gogoberidze et al. 2005). 

3. Pulsars with A/3 < They should show observed intensity variations, modulated 
with the period of the drift waves. 

4. Pulsars with A/3 > $ and (P dr — mP)/P « 1 (where m is a positive integer). These 
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appear different from the real, long observable rotation period. 

Thus, long-period radio pulsars represent one of the branches of usual pulsars and 
must be considered in the frameworks of traditional theories for the specific values of the 
parameters (Set A). 

We thank the referee for helpful suggestions. This work was partly supported by grants 
from the Georgian Academy of Science, the Russian Foundation for Basic Research (project 
03-02-16509), and the NSF (project 00-98685). 
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Ffs) 

Fig. 1. — Line A. "Death line" when curvature radiation with the sunspot configuration 
of the magnetic field is taken as the mechanism of producing the gamma quanta. Line B: 

B = Bschw 
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Fig. 2. — Geometry of rotation (fl), emission (K), and observers (A) axes. Angles 5 and i? 
are constants, while j3 and a are oscillating with time. 
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Fig. 3. — Oscillating behavior of a with time for /3 = 5 ~ 0.12, A(3 = 0.12, u^ r = 
2n/17(s- 1 ), Q = 2tt/0.85(s- 1 ), <p = 0; h = 2n/Q, t 2 = 4vr/fi 
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P(S) 



Fig. 4. — (P — B s ) diagram of the real positions of the considered pulsars, corresponding to 
Table 2. 



-18- 



Table 1: Radio pulsars with long periods. 



No 


Pulsar 


P(s) 


P(10- l5 s ■ s- 1 ) 


B S (10 12 G) 


E(10 32 erg ■ s' 1 ) 


I 


PSR J2144-3933 


8.5 


0.48 


2 


0.00032 


II 


PSR J1847-0130 


6.7 


1275 


94 


1.7 


III 


PSR J1814-1744 


4.0 


743 


55 


4.7 



Table 2. The 'real' values of pulsar parameters. 



Pulsar 


m 


Pdr 
(s) 


I' 

(s) 


P 

(10~ 15 s • s- 1 ) 


B s 
(10 12 G) 


E 

(10 32 er 9 ■ s- 1 ) 


A/3 
(deg) 


(deg) 


(deg) 


Wio/P 


PSR J2144-3933 


10 


17.0 


0.85 


0.048 


0.2 


0.032 


7° 


7° 


1.5° 


0.1 


PSR J1847-0130 


6 


13.4 


1.12 


210 


16 


61 


5° 


5° 


3° 


0.3 


PSR J1814-1744 


8 


8.0 


0.5 


93 


6.9 


300 


5° 


5° 


2° 


0.2 



